Disentangling top-down drivers of mortality underlying diel population dynamics of Prochlorococcus in the North Pacific Subtropical Gyre.

Photosynthesis fuels primary production at the base of marine food webs. Yet, in many surface ocean ecosystems, diel-driven primary production is tightly coupled to daily loss. This tight coupling raises the question: which top-down drivers predominate in maintaining persistently stable picocyanobacterial populations over longer time scales? Motivated by high-frequency surface water measurements taken in the North Pacific Subtropical Gyre (NPSG), we developed multitrophic models to investigate bottom-up and top-down mechanisms underlying the balanced control of Prochlorococcus populations. We find that incorporating photosynthetic growth with viral- and predator-induced mortality is sufficient to recapitulate daily oscillations of Prochlorococcus abundances with baseline community abundances. In doing so, we infer that grazers in this environment function as the predominant top-down factor despite high standing viral particle densities. The model-data fits also reveal the ecological relevance of light-dependent viral traits and non-canonical factors to cellular loss. Finally, we leverage sensitivity analyses to demonstrate how variation in life history traits across distinct oceanic contexts, including variation in viral adsorption and grazer clearance rates, can transform the quantitative and even qualitative importance of top-down controls in shaping Prochlorococcus population dynamics.

Widespread yet persistent low abundance of TIM5-like cyanophages in the oceans.

Ocean ecosystems are inhabited by a diverse set of viruses that impact microbial mortality and evolution. However, the distribution and abundances of specific viral lineages, particularly those from the large bank of rare viruses, remains largely unknown. Here, we assessed the diversity and abundance of the TIM5-like cyanophages. The sequencing of three new TIM5-like cyanophage genomes and environmental amplicons of a signature gene from the Red Sea revealed highly conserved gene content and sequence similarity. We adapted the polony method, a solid-phase polymerase chain reaction assay, to quantify TIM5-like cyanophages during three 2000 km expeditions in the Pacific Ocean and four annual cycles in the Red Sea. TIM5-like cyanophages were widespread, detected at all latitudes and seasons surveyed throughout the photic zone. Yet they were generally rare, ranging between <100 and 4000 viruses·ml-1 . Occasional peaks in abundance of 10- to 100-fold were observed, reaching 71,000 viruses·ml-1 . These peaks were ephemeral and seasonally variable in the Red Sea. Infection levels, quantified during one such peak, were very low. These characteristics of low diversity and abundance, as well as variable outbreaks, distinguishes the TIM5-like lineage from other major cyanophage lineages and illuminates that rare virus lineages can be persistent and widespread in the oceans.

Cyanophages from a less virulent clade dominate over their sister clade in global oceans.

Environmental virus communities are highly diverse. However, the infection physiology underlying the evolution of diverse phage lineages and their ecological consequences are largely unknown. T7-like cyanophages are abundant in nature and infect the marine unicellular cyanobacteria, Synechococcus and Prochlorococcus, important primary producers in the oceans. Viruses belonging to this genus are divided into two distinct phylogenetic clades: clade A and clade B. These viruses have narrow host-ranges with clade A phages primarily infecting Synechococcus genotypes, while clade B phages are more diverse and can infect either Synechococcus or Prochlorococcus genotypes. Here we investigated infection properties (life history traits) and environmental abundances of these two clades of T7-like cyanophages. We show that clade A cyanophages have more rapid infection dynamics, larger burst sizes and greater virulence than clade B cyanophages. However, clade B cyanophages were at least 10-fold more abundant in all seasons, and infected more cyanobacteria, than clade A cyanophages in the Red Sea. Models predicted that steady-state cyanophage abundances, infection frequency, and virus-induced mortality, peak at intermediate virulence values. Our findings indicate that differences in infection properties are reflected in virus phylogeny at the clade level. They further indicate that infection properties, together with differences in subclade diversity and host repertoire, have important ecological consequences with the less aggressive, more diverse virus clade having greater ecological impacts.

Viruses affect picocyanobacterial abundance and biogeography in the North Pacific Ocean.

The photosynthetic picocyanobacteria Prochlorococcus and Synechococcus are models for dissecting how ecological niches are defined by environmental conditions, but how interactions with bacteriophages affect picocyanobacterial biogeography in open ocean biomes has rarely been assessed. We applied single-virus and single-cell infection approaches to quantify cyanophage abundance and infected picocyanobacteria in 87 surface water samples from five transects that traversed approximately 2,200 km in the North Pacific Ocean on three cruises, with a duration of 2-4 weeks, between 2015 and 2017. We detected a 550-km-wide hotspot of cyanophages and virus-infected picocyanobacteria in the transition zone between the North Pacific Subtropical and Subpolar gyres that was present in each transect. Notably, the hotspot occurred at a consistent temperature and displayed distinct cyanophage-lineage composition on all transects. On two of these transects, the levels of infection in the hotspot were estimated to be sufficient to substantially limit the geographical range of Prochlorococcus. Coincident with the detection of high levels of virally infected picocyanobacteria, we measured an increase of 10-100-fold in the Synechococcus populations in samples that are usually dominated by Prochlorococcus. We developed a multiple regression model of cyanophages, temperature and chlorophyll concentrations that inferred that the hotspot extended across the North Pacific Ocean, creating a biological boundary between gyres, with the potential to release organic matter comparable to that of the sevenfold-larger North Pacific Subtropical Gyre. Our results highlight the probable impact of viruses on large-scale phytoplankton biogeography and biogeochemistry in distinct regions of the oceans.

Genetic engineering of marine cyanophages reveals integration but not lysogeny in T7-like cyanophages.

Marine cyanobacteria of the genera Synechococcus and Prochlorococcus are the most abundant photosynthetic organisms on earth, spanning vast regions of the oceans and contributing significantly to global primary production. Their viruses (cyanophages) greatly influence cyanobacterial ecology and evolution. Although many cyanophage genomes have been sequenced, insight into the functional role of cyanophage genes is limited by the lack of a cyanophage genetic engineering system. Here, we describe a simple, generalizable method for genetic engineering of cyanophages from multiple families, that we named REEP for REcombination, Enrichment and PCR screening. This method enables direct investigation of key cyanophage genes, and its simplicity makes it adaptable to other ecologically relevant host-virus systems. T7-like cyanophages often carry integrase genes and attachment sites, yet exhibit lytic infection dynamics. Here, using REEP, we investigated their ability to integrate and maintain a lysogenic life cycle. We found that these cyanophages integrate into the host genome and that the integrase and attachment site are required for integration. However, stable lysogens did not form. The frequency of integration was found to be low in both lab cultures and the oceans. These findings suggest that T7-like cyanophage integration is transient and is not part of a classical lysogenic cycle.

Adaptation of the polony technique to quantify Gokushovirinae, a diverse group of single-stranded DNA phage.

Advances in metagenomics have revealed the ubiquity of single-stranded DNA (ssDNA) phage belonging to the subfamily Gokushovirinae in the oceans; however, the abundance and ecological roles of this group are unknown. Here, we quantify gokushoviruses through adaptation of the polony method, in which viral template DNA is immobilized in a gel, amplified by PCR, and subsequently detected by hybridization. Primers and probes for this assay were designed based on PCR amplicon diversity of gokushovirus major capsid protein gene sequences from a depth profile in the Gulf of Aqaba, Red Sea sampled in September 2015. At ≥95% identity, these 87 gokushovirus sequences formed 14 discrete clusters with the largest clades showing distinct depth distributions. The application of the polony method enabled the first quantification of gokushoviruses in any environment. The gokushoviruses were most abundant in the upper 40 m of the stratified water column, with a subsurface peak in abundance of 1.26 × 105 viruses ml-1 . These findings suggest that discrete gokushovirus genotypes infect bacterial hosts that differentially partition in the water column. Since the designed primers and probe are conserved across marine ecosystems, this polony method can be applied broadly for the quantification of gokushoviruses throughout the global oceans.

A single-cell polony method reveals low levels of infected Prochlorococcus in oligotrophic waters despite high cyanophage abundances.

Long-term stability of picocyanobacteria in the open oceans is maintained by a balance between synchronous division and death on daily timescales. Viruses are considered a major source of microbial mortality, however, current methods to measure infection have significant methodological limitations. Here we describe a method that pairs flow-cytometric sorting with a PCR-based polony technique to simultaneously screen thousands of taxonomically resolved individual cells for intracellular virus DNA, enabling sensitive, high-throughput, and direct quantification of infection by different virus lineages. Under controlled conditions with picocyanobacteria-cyanophage models, the method detected infection throughout the lytic cycle and discriminated between varying infection levels. In North Pacific subtropical surface waters, the method revealed that only a small percentage of Prochlorococcus (0.35-1.6%) were infected, predominantly by T4-like cyanophages, and that infection oscillated 2-fold in phase with the diel cycle. This corresponds to 0.35-4.8% of Prochlorococcus mortality daily. Cyanophages were 2-4-fold more abundant than Prochlorococcus, indicating that most encounters did not result in infection and suggesting infection is mitigated via host resistance, reduced phage infectivity and inefficient adsorption. This method will enable quantification of infection for key microbial taxa across oceanic regimes and will help determine the extent that viruses shape microbial communities and ecosystem level processes.

Quantification of T4-Like and T7-Like Cyanophages Using the Polony Method Show They Are Significant Members of the Virioplankton in the North Pacific Subtropical Gyre.

The North Pacific Subtropical Gyre (NPSG) is one of the largest biomes on Earth, with the cyanobacterium Prochlorococcus being the most abundant primary producer year-round. Viruses that infect cyanobacteria (cyanophages) influence cyanobacterial mortality, diversity and evolution. Two major cyanophage families are the T4-like cyanomyoviruses and T7-like cyanopodoviruses, yet their abundances and distribution patterns remain unknown due to difficulty in quantifying their populations. To address this limitation, we previously adapted the polony method (for PCR colony) to quantify T7-like cyanophages and applied it to spring populations in the Red Sea. Here, we further adapted the method for the quantification of T4-like cyanophages and analyzed the abundances of T4-like and T7-like cyanophage populations in the photic zone of the NPSG in summer 2015 and spring 2016. Combined, the peak abundances of these two cyanophage families reached 2.8 × 106 and 1.1 × 106 cyanophages ⋅ ml-1 in the summer and spring, respectively. They constituted between 3 and 16% of total virus-like particles (VLPs), comprising a substantial component of the virioplankton in the NPSG. While both cyanophage families were highly abundant, the T4-like cyanophages were generally 1.3-4.4 fold more so. In summer, cyanophages had similar and reproducible distribution patterns with depth. Abundances were relatively low in the upper mixed layer and increased to form a pronounced subsurface peak at 100 m (1.9 × 106 and 9.1 × 105 phages ⋅ ml-1 for the T4-like and T7-like cyanophages, respectively), coincident with the maximum in Prochlorococcus populations. Less vertical structure in cyanophage abundances was apparent in the spring profile, despite a subsurface peak in Prochlorococcus numbers. In the summer upper mixed layer, cyanophages constituted a smaller proportion of VLPs than below it and cyanophage to cyanobacteria ratios were considerably lower (1.3-2.8) than those of VLPs to bacteria (8.1-21.2). Differences in abundances between the two families and their contribution to VLPs with depth suggest differences in cyanophage production and/or decay processes relative to other members of the virioplankton in the upper mixed layer. These findings highlight the importance of quantifying distinct populations within the virioplankton to gain accurate understanding of their distribution patterns.

Resistance in marine cyanobacteria differs against specialist and generalist cyanophages.

Long-term coexistence between unicellular cyanobacteria and their lytic viruses (cyanophages) in the oceans is thought to be due to the presence of sensitive cells in which cyanophages reproduce, ultimately killing the cell, while other cyanobacteria survive due to resistance to infection. Here, we investigated resistance in marine cyanobacteria from the genera Synechococcus and Prochlorococcus and compared modes of resistance against specialist and generalist cyanophages belonging to the T7-like and T4-like cyanophage families. Resistance was extracellular in most interactions against specialist cyanophages irrespective of the phage family, preventing entry into the cell. In contrast, resistance was intracellular in practically all interactions against generalist T4-like cyanophages. The stage of intracellular arrest was interaction-specific, halting at various stages of the infection cycle. Incomplete infection cycles proceeded to various degrees of phage genome transcription and translation as well as phage genome replication in numerous interactions. In a particularly intriguing case, intracellular capsid assembly was observed, but the phage genome was not packaged. The cyanobacteria survived the encounter despite late-stage infection and partial genome degradation. We hypothesize that this is tolerated due to genome polyploidy, which we found for certain strains of both Synechococcus and Prochlorococcus Our findings unveil a heavy cost of promiscuous entry of generalist phages into nonhost cells that is rarely paid by specialist phages and suggests the presence of unknown mechanisms of intracellular resistance in the marine unicellular cyanobacteria. Furthermore, these findings indicate that the range for virus-mediated horizontal gene transfer extends beyond hosts to nonhost cyanobacterial cells.

Adaptation to sub-optimal hosts is a driver of viral diversification in the ocean.

Cyanophages of the Myoviridae family include generalist viruses capable of infecting a wide range of hosts including those from different cyanobacterial genera. While the influence of phages on host evolution has been studied previously, it is not known how the infection of distinct hosts influences the evolution of cyanophage populations. Here, using an experimental evolution approach, we investigated the adaptation of multiple cyanophage populations to distinct cyanobacterial hosts. We show that when infecting an "optimal" host, whose infection is the most efficient, phage populations accumulated only a few mutations. However, when infecting "sub-optimal" hosts, different mutations spread in the phage populations, leading to rapid diversification into distinct subpopulations. Based on our results, we propose a model demonstrating how shifts in microbial abundance, which lead to infection of "sub-optimal" hosts, act as a driver for rapid diversification of viral populations.

Single cell genomes of Prochlorococcus, Synechococcus, and sympatric microbes from diverse marine environments.

Prochlorococcus and Synechococcus are the dominant primary producers in marine ecosystems and perform a significant fraction of ocean carbon fixation. These cyanobacteria interact with a diverse microbial community that coexists with them. Comparative genomics of cultivated isolates has helped address questions regarding patterns of evolution and diversity among microbes, but the fraction that can be cultivated is miniscule compared to the diversity in the wild. To further probe the diversity of these groups and extend the utility of reference sequence databases, we report a data set of single cell genomes for 489 Prochlorococcus, 50 Synechococcus, 9 extracellular virus particles, and 190 additional microorganisms from a diverse range of bacterial, archaeal, and viral groups. Many of these uncultivated single cell genomes are derived from samples obtained on GEOTRACES cruises and at well-studied oceanographic stations, each with extensive suites of physical, chemical, and biological measurements. The genomic data reported here greatly increases the number of available Prochlorococcus genomes and will facilitate studies on evolutionary biology, microbial ecology, and biological oceanography.

Quantification of diverse virus populations in the environment using the polony method.

Viruses are globally abundant and extremely diverse in their genetic make-up and in the hosts they infect. Although they influence the abundance, diversity and evolution of their hosts, current methods are inadequate for gaining a quantitative understanding of their impact on these processes. Here we report the adaptation of the solid-phase single-molecule PCR polony method for the quantification of taxonomically relevant groups of diverse viruses. Using T7-like cyanophages as our model, we found the polony method to be far superior to regular quantitative PCR methods and droplet digital PCR when degenerate primers were used to encompass the group's diversity. This method revealed that T7-like cyanophages were highly abundant in the Red Sea in spring 2013, reaching 770,000 phages ml-1, and displaying a similar depth distribution pattern to cyanobacteria. Furthermore, the abundances of two major clades within the T7-like cyanophages differed dramatically throughout the water column: clade B phages that carry the psbA photosynthesis gene and infect either Synechococcus or Prochlorococcus were at least 20-fold more abundant than clade A phages that lack psbA and infect Synechococcus hosts. Such measurements are of paramount importance for understanding virus population dynamics and the impact of viruses on different microbial taxa and for modelling viral influence on ecosystem functioning on a global scale.

A myovirus encoding both photosystem I and II proteins enhances cyclic electron flow in infected Prochlorococcus cells.

Cyanobacteria are important contributors to primary production in the open oceans. Over the past decade, various photosynthesis-related genes have been found in viruses that infect cyanobacteria (cyanophages). Although photosystem II (PSII) genes are common in both cultured cyanophages and environmental samples 1-4 , viral photosystem I (vPSI) genes have so far only been detected in environmental samples 5,6 . Here, we have used a targeted strategy to isolate a cyanophage from the tropical Pacific Ocean that carries a PSI gene cassette with seven distinct PSI genes (psaJF, C, A, B, K, E, D) as well as two PSII genes (psbA, D). This cyanophage, P-TIM68, belongs to the T4-like myoviruses, has a prolate capsid, a long contractile tail and infects Prochlorococcus sp. strain MIT9515. Phage photosynthesis genes from both photosystems are expressed during infection, and the resultant proteins are incorporated into membranes of the infected host. Moreover, photosynthetic capacity in the cell is maintained throughout the infection cycle with enhancement of cyclic electron flow around PSI. Analysis of metagenomic data from the Tara Oceans expedition 7 shows that phages carrying PSI gene cassettes are abundant in the tropical Pacific Ocean, composing up to 28% of T4-like cyanomyophages. They are also present in the tropical Indian and Atlantic Oceans. P-TIM68 populations, specifically, compose on average 22% of the PSI-gene-cassette carrying phages. Our results suggest that cyanophages carrying PSI and PSII genes are likely to maintain and even manipulate photosynthesis during infection of their Prochlorococcus hosts in the tropical oceans.

Two Synechococcus genes, Two Different Effects on Cyanophage Infection.

Synechococcus is an abundant marine cyanobacterium that significantly contributes to primary production. Lytic phages are thought to have a major impact on cyanobacterial population dynamics and evolution. Previously, an investigation of the transcriptional response of three Synechococcus strains to infection by the T4-like cyanomyovirus, Syn9, revealed that while the transcript levels of the vast majority of host genes declined soon after infection, those for some genes increased or remained stable. In order to assess the role of two such host-response genes during infection, we inactivated them in Synechococcus sp. strain WH8102. One gene, SYNW1659, encodes a domain of unknown function (DUF3387) that is associated with restriction enzymes. The second gene, SYNW1946, encodes a PIN-PhoH protein, of which the PIN domain is common in bacterial toxin-antitoxin systems. Neither of the inactivation mutations impacted host growth or the length of the Syn9 lytic cycle. However, the DUF3387 mutant supported significantly lower phage DNA replication and yield of phage progeny than the wild-type, suggesting that the product of this host gene aids phage production. The PIN-PhoH mutant, on the other hand, allowed for significantly higher Syn9 genomic DNA replication and progeny production, suggesting that this host gene plays a role in restraining the infection process. Our findings indicate that host-response genes play a functional role during infection and suggest that some function in an attempt at defense against the phage, while others are exploited by the phage for improved infection.

Genetic hurdles limit the arms race between Prochlorococcus and the T7-like podoviruses infecting them.

Phages and hosts coexist in nature with a high degree of population diversity. This is often explained through coevolutionary models, such as the arms race or density-dependent fluctuating selection, which differ in assumptions regarding the emergence of phage mutants that overcome host resistance. Previously, resistance in the abundant marine cyanobacterium, Prochlorococcus, was found to occur frequently. However, little is known about the ability of phages to overcome this resistance. Here we report that, in some cases, T7-like cyanophage mutants emerge to infect resistant Prochlorococcus strains. These resistance-breaking phages retained the ability to infect the wild-type host. However, fitness of the mutant phages differed on the two hosts. Furthermore, in one case, resistance-breaking was accompanied by costs of decreased fitness on the wild-type host and decreased adsorption specificity, relative to the wild-type phage. In two other cases, fitness on the wild-type host increased. Whole-genome sequencing revealed mutations in probable tail-related genes. These were highly diverse in isolates and natural populations of T7-like cyanophages, suggesting that antagonistic coevolution enhances phage genome diversity. Intriguingly, most interactions did not yield resistance-breaking phages. Thus, resistance mutations raise genetic barriers to continuous arms race cycles and are indicative of an inherent asymmetry in coevolutionary capacity, with hosts having the advantage. Nevertheless, phages coexist with hosts, which we propose relies on combined, parallel action of a limited arms race, fluctuating selection and passive host-switching within diverse communities. Together, these processes generate a constantly changing network of interactions, enabling stable coexistence between hosts and phages in nature.

A Novel Strategy for Exploitation of Host RNase E Activity by a Marine Cyanophage.

Previous studies have shown that infection of Prochlorococcus MED4 by the cyanophage P-SSP7 leads to increased transcript levels of host endoribonuclease (RNase) E. However, it has remained enigmatic whether this is part of a host defense mechanism to degrade phage messenger RNA (mRNA) or whether this single-strand RNA-specific RNase is utilized by the phage. Here we describe a hitherto unknown means through which this cyanophage increases expression of RNase E during phage infection and concomitantly protects its own RNA from degradation. We identified two functionally different RNase E mRNA variants, one of which is significantly induced during phage infection. This transcript lacks the 5' UTR, is considerably more stable than the other transcript, and is likely responsible for increased RNase E protein levels during infection. Furthermore, selective enrichment and in vivo analysis of double-stranded RNA (dsRNA) during infection revealed that phage antisense RNAs (asRNAs) sequester complementary mRNAs to form dsRNAs, such that the phage protein-coding transcriptome is nearly completely covered by asRNAs. In contrast, the host protein-coding transcriptome is only partially covered by asRNAs. These data suggest that P-SSP7 orchestrates degradation of host RNA by increasing RNase E expression while masking its own transcriptome from RNase E degradation in dsRNA complexes. We propose that this combination of strategies contributes significantly to phage progeny production.

Variability in progeny production and virulence of cyanophages determined at the single-cell level.

Little information regarding viral progeny production (burst size) and host mortality (viral virulence) is currently available for environmentally relevant phages. This is partially due to the difficulty in accurately measuring these infection properties with existing methods. Here, we set up a simple system for determining viral virulence and burst size at the single-cell level following flow cytometric separation of infected cells. We applied this assay to two distinct cyanomyoviruses, Syn9 and S-TIM5, during infection of two marine Synechococcus strains each. We found that virulence ranged from 44%-82%, differing for the same phage on different hosts. Average burst sizes ranged from 21-43 infective viruses/cell, and differed with host for Syn9, whereas the burst size of S-TIM5 was similar for both hosts. In addition, virulence and burst sizes were different for the two phages when infecting their common host. Furthermore, wide-ranging cell-to-cell variability was found for single-cell burst sizes in each of the four interactions, ranging from 2 to over 100 infective viruses/cell. This variability, discerned at both the population and single-cell levels under controlled laboratory conditions, is likely to be much more complex in natural environments.

Transcriptome dynamics of a broad host-range cyanophage and its hosts.

Cyanobacteria are highly abundant in the oceans and are constantly exposed to lytic viruses. The T4-like cyanomyoviruses are abundant in the marine environment and have broad host-ranges relative to other cyanophages. It is currently unknown whether broad host-range phages specifically tailor their infection program for each host, or employ the same program irrespective of the host infected. Also unknown is how different hosts respond to infection by the same phage. Here we used microarray and RNA-seq analyses to investigate the interaction between the Syn9 T4-like cyanophage and three phylogenetically, ecologically and genomically distinct marine Synechococcus strains: WH7803, WH8102 and WH8109. Strikingly, Syn9 led a nearly identical infection and transcriptional program in all three hosts. Different to previous assumptions for T4-like cyanophages, three temporally regulated gene expression classes were observed. Furthermore, a novel regulatory element controlled early-gene transcription, and host-like promoters drove middle gene transcription, different to the regulatory paradigm for T4. Similar results were found for the P-TIM40 phage during infection of Prochlorococcus NATL2A. Moreover, genomic and metagenomic analyses indicate that these regulatory elements are abundant and conserved among T4-like cyanophages. In contrast to the near-identical transcriptional program employed by Syn9, host responses to infection involved host-specific genes primarily located in hypervariable genomic islands, substantiating islands as a major axis of phage-cyanobacteria interactions. Our findings suggest that the ability of broad host-range phages to infect multiple hosts is more likely dependent on the effectiveness of host defense strategies than on differential tailoring of the infection process by the phage.

Convergent evolution toward an improved growth rate and a reduced resistance range in Prochlorococcus strains resistant to phage.

Prochlorococcus is an abundant marine cyanobacterium that grows rapidly in the environment and contributes significantly to global primary production. This cyanobacterium coexists with many cyanophages in the oceans, likely aided by resistance to numerous co-occurring phages. Spontaneous resistance occurs frequently in Prochlorococcus and is often accompanied by a pleiotropic fitness cost manifested as either a reduced growth rate or enhanced infection by other phages. Here, we assessed the fate of a number of phage-resistant Prochlorococcus strains, focusing on those with a high fitness cost. We found that phage-resistant strains continued evolving toward an improved growth rate and a narrower resistance range, resulting in lineages with phenotypes intermediate between those of ancestral susceptible wild-type and initial resistant substrains. Changes in growth rate and resistance range often occurred in independent events, leading to a decoupling of the selection pressures acting on these phenotypes. These changes were largely the result of additional, compensatory mutations in noncore genes located in genomic islands, although genetic reversions were also observed. Additionally, a mutator strain was identified. The similarity of the evolutionary pathway followed by multiple independent resistant cultures and clones suggests they undergo a predictable evolutionary pathway. This process serves to increase both genetic diversity and infection permutations in Prochlorococcus populations, further augmenting the complexity of the interaction network between Prochlorococcus and its phages in nature. Last, our findings provide an explanation for the apparent paradox of a multitude of resistant Prochlorococcus cells in nature that are growing close to their maximal intrinsic growth rates.

Host-dependent differences in abundance, composition and host range of cyanophages from the Red Sea.

Cyanobacteria coexist in the oceans with a wealth of phages that infect them. While numerous studies have investigated Synechococcus phages, much less data are available for Prochlorococcus phages. Furthermore, little is known about cyanophage composition. Here, we examined the abundance and relative composition of cyanophages on six cyanobacterial hosts in samples collected during spring and summer from the Red Sea. Maximal abundances found on Synechococcus of 35 000 phages/ml are within ranges found previously, whereas the 24 000 phages/ml found on Prochlorococcus are approximately 10-fold higher than previous findings. T7-like, T4-like and 'unknown' phages were isolated on all hosts, including many T4-like phages on high-light adapted Prochlorococcus strains, whereas TIM5-like phages were found only on Synechococcus. Large differences in cyanophage abundance and composition were found for different hosts on the same sampling date, as well as for the same host on different dates, with few predictable patterns discerned. Host range analyses showed that T7-like and TIM5-like phages were quite host-specific, whereas the breadth of hosts for T4-like phages was related to host type: those isolated on high-light adapted Prochlorococcus were considerably more host-specific than those on low-light adapted Prochlorococcus or Synechococcus. These host-related differences likely contribute to the complexity of host-phage interactions in the oceans.